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Abstract 

We present results from Chandra and XMM-Newton observations of the bright 
group of galaxies HCG 62. There are two cavities at about 30" northeast and 20" 
southwest of the central galaxy in the Chandra image. The energy spectrum shows no 
significant change in the cavity compared with that in the surrounding region. The 
radial X-ray profile is described by a sum of 3-f3 components with core radii about 
2, 10, and 160 kpc, respectively. We studied radial distributions of temperature and 
metal abundance with joint spectral fit for the Chandra and XMM-Newton data, and 
two temperatures were required in the inner r < 2 / (35 kpc) region. The sharp drop 
of temperature at r ~ 5 / implies the gravitational mass density even lower than the 
gas density, suggesting the gas may not be in hydrostatic equilibrium. Fe and Si 
abundances are 1-2 solar at the center and drop to about 0.1 solar at r ^ 10 / . O 
abundance is less than 0.5 solar and shows a flatter profile. Observed metal distri- 
bution supports the view that iron and silicon are produced by type la supernova in 
the central galaxy, while galactic winds by type II supernova have caused wide distri- 
bution of oxygen. The supporting mechanism of the cavity is discussed. Pressure for 
the sum of electrons and magnetic field is too low to displace the hot group gas, and 
the required pressure due to high energy protons are nearly 700 times higher than 
the electron pressure. This leaves the origin of the cavities a puzzle, and we discuss 
other possible origins of the cavities. 
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1 



1. Introduction 



Groups of galaxies hold significantly less amount of hot gas compared with rich clusters 
on the average, and the apparent deficiency of baryons in these low-mass systems is a problem 
in explaining the baryon budget in rich systems in terms of hierarchical merging scenario 
(Voit 2005). The observations of compact groups of galaxies are important to search for the 
hidden form of baryons and their release mechanism, which would be strongly connected with 
dynamics of gas and galaxies. In particular, some groups dominated by bright central galaxies 
have shown gas features which are strongly affected by the activity of the central galaxies. High 
sensitivity X-ray observations of groups of galaxies are powerful method to look into the role 
of central galaxies in terms of metal distribution and gas morphologies. 

The striking gas features most likely caused by the central galaxies are the X-ray cavities. 
The cavities are circular regions showing a significant depression of X-ray surface brightness. 
Nearly 20 cavities have been recognized in clusters and groups with high resolution images 
taken by ROSAT and Chandra (Birzan et al. 2004, hereafter B04; Dunn & Fabian 2004; Dunn 
et al. 2005, hereafter D05). They are located typically at 10-30 kpc from the central galaxies, 
and strong correlation with radio lobes are seen in about 10 systems. Remarkable cases are seen 
in the Perseus (e.g., Boehringer et al. 1993; Fabian et al. 2000) and Hydra A (e.g. McNamara et 
al. 2000) clusters, both showing strong correlation with the 1.4 GHz radio lobes. Several giant 
elliptical galaxies with radio robes, e.g., M84 (Finoguenov & Jones 2001), NGC 4636 (Ohto 
at al. 2003), are also known as nesting X-ray cavities. The remaining half of the cavities, on 
the other hand, are not associated with radio lobes, and they are designated as ghost cavities. 
The one in A 2597 (McNamara et al. 2001) or the outer depressions in Perseus (Fabian et 
al. 2000) are examples. Cavities are thought to be produced by jets or buoyant bubbles which 
are directly connected with the activity of central radio galaxies. 

A subsonic displacement of the gas would create a low density, rising bubble keeping the 
pressure balance with the surrounding ICM. It appears to be in general supposed for cavities 
that non-thermal pressure originated in relativistic particles and/or magnetic fields in the radio 
lobe is probably large enough to balance with the surrounding ICM gas pressure (Fabian et 
al. 2002). This pseudo-pressure balance is justified by the fact that there are no evidence for 
shock-heated gas around the radio lobes in almost all of the X-ray cavities observed so far, except 
for MKW 3s (Mazzotta et al. 2002). This general scenario has been modeled theoretically, and 
has at least qualitatively reproduced the morphology of cavities (e.g., Churazov et al. 2001). 
Energy density of relativistic electrons inferred from the synchrotron radio emission is almost 
always smaller than that required to offset the hot gas by orders of magnitude, and it is discussed 
that energy density of protons are higher than those due to electrons by factors of 100-1000 
(D05). However, there is no direct evidence indicating that such a high energy density is really 
carried by protons. This situation is the severest in the case of ghost cavities. In this view, it 
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is important to examine ghost cavities in groups of galaxies where the gas is relatively cool and 
non-thermal effect can be recognized somewhat easily. 

Detailed studies on the metal distribution in clusters and groups have been carried out 
using ASCA, BeppoSAX, Chandra and XMM-Newton. Distribution of iron and silicon indicate 
strong central concentration in clusters and groups characterized by bright central galaxies, and 
the excess iron mass is found to correlate with the luminosity of the cD galaxy (De Grandi et 
al. 2004). This indicates that iron and silicon (main products from type la supernova; SN la) 
trace the enhanced star-formation activity in bright galaxies. On the other hand, distribution of 
oxygen (i.e. type II supernova product; SN II) is not well understood. Matsushita et al. (2003) 
showed that oxygen distribution around M 87 is flatter than those of iron and silicon, with the 
level about half as much as the others. Such low oxygen abundances are also derived in the 
centers of other clusters and groups (e.g. Buote et al. 2003b; Xue et al. 2004). For the study of 
oxygen distribution, low temperature systems such as groups of galaxies are suitable targets. 

In this paper, we report the results from Chandra and XMM-Newton observations of 
HCG 62, which is one of the nearest Hickson compact galaxy groups (Hickson et al. 1989). The 
whole group consists of 63 galaxies (Mulchaey et al. 2003) within a radius of 50 / (900 kpc), but 
the central region is dominated by 4 galaxies. HCG 62 is the brightest group of galaxies in the 
X-ray band, and the extended X-ray emission was first discovered by Ponman & Bertram (1993) 
from the ROSAT PSPC observation. Based on the ASCA observation, Fukazawa et al. (2001) 
detected excess hard X-ray emission, and Finoguenov & Ponman (1999) report strong central 
concentration of iron. Using the high resolution image of Chandra, Vrtilek et al. (2001,2002) 
detected two ghost cavities, which is the first report of cavities in groups of galaxies. 

This paper is organized as follows: In § 2 we describe the Chandra and XMM-Newton 
observations and the data reduction. In § 3 we give the image of HCG 62 of both Chandra 
and XMM-Newton, in §4 we describe the X-ray cavity structure using Chandra image. In 
§ 5 we present our results on the temperature profiles and the abundances profiles of Fe and 
a-elements (Si, Mg, and O). In §6-8 we give discussions of the obtained results, and finally we 
summarize our conclusions in §9. Throughout this paper we adopt = 1 _ = 0.73 and 
/i 70 = H /{70 km s _1 Mpc _1 ) = 1; V corresponds to 17.8 kpc at z — 0.0145. The quoted errors 
indicate the 90% confidence range, unless otherwise stated. We use the solar abundance ratio 
of Anders and Grevesse (1989). 

2. Observation and Data Reduction 

2.1. Chandra Observation 

HCG 62 was observed on 25 January 2000 with the Advanced CCD imaging Spectrometer 
(ACIS) 12, 13, S2, S3, and S4 chips operated at the CCD temperature of — 110°C with a frame 
readout time of 3.24 s. We used the data of only the ACIS-S3 chip (8.4 / x 8.4') covering the cen- 
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tral part of HCG 62 in this paper. The pointing coordinates were (12 h 59 m 05?70, -09°12'20"00) 
(J2000) and the total exposure was 49.15 ks. According to the "Chandra Aspect" web page, 
http://cxc.harvard.edu/cal/ASPECT/, the astrometry offset of RA = -0.03", Dec = 0.22" has 
been corrected for in the data, providing a radius of 0.6 arcsec as the absolute position accu- 
racy in 90% confidence. The CCD temperature of ACIS was reduced to — 120°C soon after the 
observation of HCG 62 in January 2000 due to the increase of charge transfer inefficiency (CTI) 
caused by the radiation damage in orbit (Grant et al. 2005). The data were telemetered in the 
Faint mode, and events with the ASCA grades of 0, 2, 3, 4, and 6 were used. Bad pixels, bad 
columns, and the columns next to bad columns and to the chip node boundaries are excluded. 
In order to remove periods of anomalous background levels, we further filtered the events using 
the 0.3-10 keV band light curve of the whole ACIS-S3 chip in 200 s bin, and discarded periods 
which exceeded by 3a above the mean quiescent rate of 4.5 c/s/chip. The net exposure time 
after the screening was 48,013 s. The data reduction was performed using CIAO version 3.1 
with CALDB version 2.29, and the spectral fitting was done by XSPEC version 11.3.0t. 

All of the X-ray spectra were extracted using the pulse-height invariant (PI) values, 
which were recomputed using the latest gain file acisD1999-09-16gainN0005 . f its, appropri- 
ate for the CCD temperature of — 110°C. In the spectral fitting, we initially generated the 
response matrix file (RMF) of the ACIS-S3 using the CIAO "mkrmP task with the input FEF 
(FITS Embedded Function) file of acisD1999-09-16f ef _phaN0002 . f its, which was chosen by 
default. However, we found that the Si line of the IGM emission was significantly broader than 
the response (a = 42t§ eV), which caused the fit statistics not acceptable. In order to examine 
whether this broad Si line is a target specific issue or a common calibration problem, we fur- 
ther checked the neutral Si line, K al = 1739.98 eV and K a2 = 1739.38 eV, which is originated 
in the instrumental background of the CCD, using the blank-sky data obtained at — 110°C 
(figure 1 (a)). The fit result gave the Gaussian a = 29 ±4 eV (FWHM = 68 ±9 eV) and the 
line center energy of 1753 ±3 eV, which is significantly broad, too. We splitted the blank-sky 
spectrum into two, dividing the integration region into half, i.e. upper rows and lower rows 
of the CCD, although the results were similar. This result indicates that the charge transfer 
inefficiency (CTI) correction is working well. 

We therefore concluded that the broad Si line was resulted because the generated RMF 
had too sharp Gaussian core. In fact, the derived Si line width for the blank-sky is consistent 
with a calibration document, ACIS Memo #182 by LaMarr (2000). The FWHM of the Gaussian 
core for the RMF denoted in the FEF file is G1_FWHM = 73 eV for ACIS-S3 (CCD_ID=7) at 
the energy of 1.8 keV, while it should be around 125 eV according to the memo. Our blank-sky 
data indicated \/73 2 + 68 2 = 100 eV, which is slightly better than the value in the LaMarr's 
memo, possibly due to the updated gain file. We therefore created a new FEF file to adopt 
broader Gaussian core over the whole energy range, modifying the G1_FWHM column into 
(16 x ENERGY + G1_FWHM). Using the new FEF file, we generated the RMF files with the 
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"mkrmf" task for the following analysis. 

Regarding the "blank-sky" background, we adopted a series of observations compiled 
by Markevitch (http:/ /cxc. harvard.edu/contrib/maxim/acisbg/), when the CCD temperature 
was — 110°C. The "blank-sky" background events were also screened to remove background 
flares in the same way as described above. The additional component due to the Galactic soft 
background are considered in §2.3. 

2.2. XMM-Newton observation 

The XMM-Newton observation of HCG 62 was carried out on 15 January 2003, assigned 
for 12.6 ks. The EPIC cameras were operated in full-frame mode for MOS, and in extended- 
full-frame mode for pn. The medium filters were used for both cameras. XMM-Newton covers 
the wider field of view of r < 15 / than Chandra, while the angular resolution of 15" half-power 
diameter is broader. Data reduction was performed using SAS version 6.0, and the spectral 
fitting was done by XSPEC version 11.3.0t. We selected events with pattern 0-12 for MOS and 
0-4 for pn, and flag = for both. Bad pixels and bad columns were excluded. We calculated 
the count rate distribution with 100 s intervals over the 0.3-10 keV range using the whole chip 
of each sensor, and rejected periods by requiring all the count rates to be within ±2<r around 
the mean. We iterated the process until the number of rejected data in a step reached less than 
5% of the 2a compared with the previous value. Because we did not find large time variation 
nor flares, the net exposure time after the screening was 12.6, 12.5, and 9.2 ks for MOS1, MOS2, 
and pn, respectively. 

To correct for the vignetting effect, the SAS "evigweight" task was applied for each 
event file. All the X-ray spectra were extracted using the vignetting-weighted events. The 
response files for spectral fittings were generated in the standard way with the SAS "rmfgen" 
and "arfgen" tasks at on-axes of the X-ray mirrors. The MOS1 and MOS2 spectra are summed 
up, and the energy range of 1.4367-1.5367 keV, where the background Al-K a (1.4867 keV) is 
strong as seen in figure 1 (b) and (c), is ignored in the spectral fitting for both MOS and pn. 

The background event dataset created by Read & Ponman (2003) was adopted as the 
"blank-sky" background for the XMM-Newton observation. We applied the same selection 
criteria as the source events to the "blank-sky" background. The "blank-sky" background was 
scaled to 0.94 for MOS1, 1.00 for MOS2, and 0.84 for pn, respectively. These factors were 
determined by the count rate ratio within r < 14' from the central galaxy HCG 62a in 10- 
12 keV (MOS) or 12-14 keV (pn) band, where the instrumental background is dominant for 
each sensor. Katayama et al. (2004) and Nevalainen et al. (2005) have studied the background 
data of XMM-Newton in detail, and it is reported that the the 90% confidence background 
uncertainty is ±5% in 4-7 keV and ±20% in 0.8-1 keV. It is confirmed that our results does 
not change significantly within errors when this level of uncertainty is considered. 
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(a) ACIS-S3 (b) M0S1+2 (c) pn 




Energy (keV) Energy (keV) Energy (keV) 

Fig. 1. X-ray spectra of (a) Chandra ACIS-S3, (b) XMM-Newton MOS1+2, and (c) pn, in 0.5-4 keV 
within r < 2' from the central galaxy HCG 62a, subtracted by the "blank-sky" background. The 
"blank-sky" spectra are indicated by gray lines. The additional background components of the 0.3 keV 
Mekal and the T = 1.5 power-law models are shown by dashed black lines, and their sum is indicated by 
solid black lines. 

2.3. Background estimation 

HCG 62 is located near the edge of the North Polar Spur (NPS), which is a large soft 
X-ray Galactic structure. The ROSAT All-Sky Survey (RASS) image at 3/4 keV band in this 
region indicates a soft X-ray excess around the position of HCG 62, although significant fraction 
of the excess is probably caused by HCG 62 itself. This emission looks more extended than 
the Chandra and XMM-Newton fields of view in the previous observations with ROSAT and 
ASCA. Therefore, we cannot use the events in outer region as the background, and we need to 
estimate the influence of the soft X-ray background to our HCG 62 data using the "blank-sky" 
background. 

There are four fields near HCG 62 already observed with Chandra, Q1246— 0542, 
NGC 4697, NGC 4594, and NGC 4782. We found that the intensities below 2 keV in source-free 
regions in these fields are higher than the level of the "blank-sky" background, while the hard 
band intensities show no excess. By fitting the energy spectra of the soft excess component with 
the XSPEC Mekal model, we obtained the temperature to be ~ 0.3 keV, consistent with the pre- 
vious result for the Galactic soft emission from NPS (Inoue et al. 1980). The surface brightness 
in the four nearby fields differ significantly, ranging in (5 10) x 10 -10 photons cm -2 s _1 arcsec -2 
in 0.5-1 keV band, and we need to estimate the soft X-ray background component based on 
the data of the HCG 62 field. 

We examined the radial profile at r < 15 / of the XMM-Newton data in 0.5-1 keV band by 
fitting with a sum of a double (3 model and a constant intensity representing the soft X-ray back- 
ground. In fitting the profile, parameters for the double (3 component were fixed at the ROSAT 
values (Zabludoff & Mulchaey 2000), and only the level of the constant component was varied 
as a free parameter. The resultant surface brightness is 2 x 10 -11 photons cm -2 s _1 arcsec -2 
in 0.5-1 keV band. In the following analysis, the excess soft component is assumed to be a 
thermal emission with the temperature of 0.3 keV and abundance of 1 solar, with the normal- 
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ization fixed at the level of the radial profile fit above. We do not apply interstellar absorption 
for this component, since its origin is supposed to be a nearby region in our galaxy. When 
the normalization of the 0.3 keV component is varied in the spectral fit, it agrees with the 
level derived from the radial profile, within the 90% confidence limit for both Chandra and 
XMM-Newton. In fact, the estimated additional soft component is much fainter than those in 
the four nearby fields, even less than the "blank-sky" background component. The "blank-sky" 
background also contains a certain fraction of the Galactic soft emission, and the HCG 62 field 
appears to require only the same level of soft X-ray background. We have also confirmed that 
our result did not change significantly within errors even doubling the intensity of the estimated 
soft background component. 

As for the additional hard X-ray background component, we take into account the excess 
emission detected with ASCA (Fukazawa et al. 2001). In the 2-10 keV range, the observed flux 
of the r = 1.5 power-law component with ASCA is 1.0 x 10 -12 ergs cm -2 s _1 at the ring-like 
area with inner and outer radii to be 5 / and 15 / . We assumed the same power-law spectrum 
in our analysis, with the spatial distribution uniform over the whole HCG 62 field for both 
Chandra and XMM-Newton. 

In the following analysis, the soft and the hard components were added as additional 
spectral models with fixed parameters. We however confirmed that the results did not vary 
significantly at the 90% confidence level, even if we did not apply these backgrounds. In figure 1, 
the soft and the hard X-ray background components as well as the "blank sky" background are 
compared with the IGM spectrum for each instrument. 

3. X-ray Image 

3.1. Chandra image 

The Chandra image of HCG 62 taken with ACIS-S3 in the 0.5-4 keV energy band is 
shown in figure 2 (a). We corrected for background and exposure, and applied Gaussian smooth- 
ing. The positional dependence of the telescope and the detector responses were corrected with 
an exposure map. The IGM emission is clearly observed, and it is extended around the central 
galaxy HCG 62a. The brightest region has a radius of about V from HCG 62a. 

We searched for discrete X-ray sources in the ACIS-S3 field using the CIAO "wavdetect" 
task. Choosing a significance parameter of 10 -6 for images in the energy bands 0.5-10 keV, 
0.5-2.0 keV, and 2.0-10.0 keV, we detected 50 sources in total, including HCG 62a, 62b, 62c, 
and 62-66. We have detected all the point sources previously identified in the SEXSI catalog 
(Harrison et al. 2003) using the same Chandra data. Optical properties of the member galaxies 
within the Chandra ACIS-S3 field of view are summarized in table 1, and the optical image is 
shown in figure 2 (b). The X-ray detected positions and the ACIS-S3 count are summarized in 
table 2. The X-ray position of HCG 62a and HCG 62b well agree with the optical coordinate 
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Table 1. Optical properties of HCG 62 and member galaxies. 



Object Optical coords. (J2000) * z * Diameter + B + Type * 

R.A. Dec. (arcsec) (mag) 

HCG 62 12 h 53 m 06?l -09°12'16"3 0.0145 Group 

HCG 62a 12 h 53 m 05f6 -09°12 / 13 // 0.0143 59.9 13.79 E3 / SO 

HCG 62b 12 h 53 m 04?4 -09°ll / 59 // 0.0119 42.9 14.21 SO / SO 

HCG 62c 12 h 53 m 09?7 -09°11 / 51 // 0.0148 40.2 15.00 SO / E 

HCG 62d 12 h 53 m 06?6 -09°15'26" 0.0136 19.7 16.30 E2 / - 

HCG 62-66 § 12 h 52 m 48?7 -09°13 / 22 // 0.0161 — / - 



Optical coordinates and redshift z by Mulchaey et al. (2003) for HCG 62 and by Zabludoff & Mulchaey (2000) for member galaxies, 
t B band effective diameter, D B =^/A B Jtt, and magnitude within \±b = 24.5 mag arcsec -2 isophote by Hickson et al. (1989). 
* Hubble morphological type classification by Hickson et al. (1989) or Shimada et al. (2000). 
§ Because HCG 62-66 is not catalogued by Hickson et al. (1989), Diameter, B, and Type are left blank. 



Table 2. X-ray properties of detected galaxies by Chandra ACIS-S3. 



Object X-ray position (J2000) * Pos. cliff, t Obs. count * Extended 

R.A. Dec. (arcsec) (cts) Y/N 

HCG 62a 12 h 53 m 05.63 -09°12 / 13. // 7 0.8 4158 ±65 Y 

HCG 62b 12 h 53 m 04.43 -09°11'59."4 0.6 1344 ±37 Y 

HCG 62c 12 h 53 m 09?76 -09°11'55."7 4.8 331 ± 19 Y? 

HCG 62d § 36 ± 7 Y ? 

HCG 62-66 12 h 52 m 48f57 -09°13'25."7 4.2 30 ±7 Y? 



* Detected position by the CI AO "wavdetect" task, with 90% confidence position accuracy of 0.6" radius, 
t Positional difference between X-ray and optical. 

* Observed count within a radius of 10 /x in 0.5-4 keV including IGM emission but subtracted by the "blank-sky". 
§ HCG 62d was not detected by the "wavdetect" task, so that the optical position was utilized. 
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Fig. 2. (a) Gaussian smoothed Chandra ACIS-S3 X-ray image in the 0.5-4 keV band. The smoothing 
scale is a = 1.5" and the image is corrected for exposure and background, (b) STScI Digitized Sky Survey 
(DSS) image of HCG 62. The ACIS-S3 field is indicated by a square, and the member galaxies are denoted 
by a, b, c, d, and 66. Coordinates of images are J2000. 

given in Zabludoff k Mulchaey (2000) within < 1". However, HCG 62c and HCG 62-66 indicate 
different positions by about 5" and 4". This is due probably to the extent of the galaxies and/or 
to the gas stripping, because the probability of miss-identification of background or foreground 
object within 5" radius is calculated to be only 1.5%. 

The X-ray emission of HCG 62a is dominated by the extended gas, for which we will 
describe the properties later. The HCG 62b galaxy was also found to be significantly more 
extended than the point spread function of Chandra. We have fitted the radial brightness 
profile, 5(r), of HCG 62b with a single (3 model, 

-3/3+1/2 



S(r) = S l + (r/R c y 



1) 



and obtained a core radius, R c = 1.5t l'l arcsec, (3 = O.5^o!i 5 an d the normalization, So = 0.8toi x 
10 -6 photons cm -2 s _1 arcsec -2 . There is also an indication for the HCG 62c galaxy to be 
extended in the Chandra image in the direction of the optical major axis (north - south), 
although it is not clear whether it is due to the extended IGM emission because of the poor 
statistics of HCG 62c. The HCG 62d galaxy was not detected by the "wavdetect" task, however, 
the "blank-sky" subtracted count exhibited a clear excess when integrated within 10" around 
the optical position as shown in table 2. It is supposed that the extended feature of HCG 62d 
may have hampered the detection by "wavdetect" . 

In the following analysis, the detected sources were excluded except for HCG 62a and 
HCG 62d. Since HCG 62a and HCG 62b are closely located, we need to separate the HCG 62b 
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10" 8 10~ 7 10" 6 10 5 count cm" 2 s" 1 arcsec 2 "10 " 1 ° 10" 9 10" 8 10 " 7 count cm" 2 s" 1 arcsec 2 



Fig. 3. (a) Same as figure 2, but the central region of HCG 62 is expanded. The detected point sources 
other than HCG 62a, 62b and 62c are excluded. The central positions for the 2-dimensional 2-/3 fit are 
denoted by cross and plus corresponding to the narrower and the wider /3 components. The red circles 
indicate a radius of V from HCG 62a and 10" from HCG 62b, respectively. The overlaid contour represent 
a radio intensity map at 1.4 GHz with 45" FWHM resolution, derived from the NRAO VLA Sky Survey 
(NVSS; Condon et al. 1998). (b) Combined image with MOS1 and MOS2 detectors in the 0.5-4 keV band 
(J2000). The image is Gaussian smoothed with a = 10" and corrected for exposure and background. A 
white square represents the ACIS-S3 field. 

component in order to examine the IGM structure. It is notable that the redshift of HCG 62b 
is smaller than the group redshift of z = 0.0145 by 0.0026 (~ 11 Mpc), hence HCG 62b is 
probably not interacting with the IGM around the group core. We produced two radial profiles 
in 0.5-4 keV band centered on HCG 62a and HCG 62b, and fitted them with a double f3 
(HCG 62a) and a single (3 (HCG 62b) models, respectively. With these fits, we evaluated the 
relative intensity of HCG 62b compared with the IGM emission as a function of the distance 
from HCG 62b. The contamination from HCG 62b becomes less than 10% of the IGM intensity 
when one goes to outside of a radius of 10". Thus, we extracted a circle with 10" radius around 
HCG 62b and carried out the analysis for the IGM emission. 

The central HCG 62 image is shown in figure 3 (a). This image indicates surface bright- 
ness depressions in the northeast and southwest directions at 20 // -30 // from HCG 62a, the 
so-called cavities as reported by Vrtilek et al. (2001,2002). We look into the properties of these 
cavities in 8 4. 
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Table 3. Best fit parameters of the radial surface brightness profiles with ACIS-S3 (r < 4') and MOS1 (1/ < r < 14' ) in 
0.5-4 keV by the 3- (3 model. First and second components were constrained to have common (3. The parameters of (3 and 
R c of the third component were fixed at the ROSAT result. 



1 2 3 



So (1(T 6 photons cm- 2 s" 1 arosec" 2 ) 0.83±g;Jg 0.29±0.03 0.0018±0.0001 

P 0.65±0.02 <- (fixed) 0.63 (fixed) 

Rc (arcmin/kpc) 0.10±g;g?/1.7±g-j 0.48t£;^/8.5 ± 0.6 9.00/159.08 (fixed) 

X 2 /dof 755.39/570 



Table 4. Best fit parameters by the 2-dimensional 2- (3 model for the central region (r < V) of HCG 62 with ACIS-S3 in 
0.5-4 keV. 



1 (narrower) 



2 (wider) 



So (10" 6 photons cm -2 s" 1 arcsec" 2 ) 0.79to:?6 

P 0.87™ 

Rc (arcmin/kpc) 0.07±g;g|/1.3±g;§ 

R.A. (J2000) 12 h 53 m 05?63 ±0.4" 

Dec. (J2000) -09°12 / 15. // 5 ±0.4" 

X 2 /dof 132.07/122 



0.63±0.06 
0.44±0.01 
0.18±g;{£/3.1±0.3 
12 h 53 m 05?90 ±0.4" 
-09°12 / ll. // 6 ±0.3" 



3.2. XMM-Newton image 

A combined X-ray image taken with MOS1 and MOS2 detectors of XMM-Newton in the 
0.5-4 keV energy band is shown in figure 3 (b). We corrected for background and exposure, and 
performed Gaussian smoothing to the image. The IGM emission is very extended beyond the 
boundary of the Chandra ACIS-S3 chip shown by a white square. The cavities are not clearly 
seen because of the poorer spatial resolution. Many point sources were detected in the outer 
region by the SAS "edetect" task, but none of them outside of r = 100" were member galaxies 
of this compact group. Point sources inside of 100" from HCG 62a were excluded using the 
Chandra data (§3.1), and the XMM-Newton data were used to mask out the sources outside 
of 100". 

3.3. Surface brightness profile 

We plot the radial surface brightness profile for the ACIS-S3 (r < 4') and MOS1 (V < 
r < 14') around HCG 62a in the energy range of 0.5-4 keV in figure 4. In this analysis, only 
the MOS1 data were used for the XMM-Newton observation because of its low background 
and good spatial resolution. HCG 62b was excluded with a circle of 10" radius. Background 
subtraction was carried out separately for each instrument. Mulchaey & Zabludoff (1998) report 
an acceptable fit with a double (3 model for the ROSAT data. As the first attempt, we fitted 
the radial profile with the same ROSAT model, however, a large discrepancy was found in the 
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Fig. 4. In the upper panel, radial profiles of the surface brightness of HCG 62 in 0.5-4 keV around the 
central galaxy HCG 62a are plotted for ACIS-S3 (r < 4') and MOS1 (1' < r < 14'), as indicated by grey and 
black lines, respectively. The orange solid line represents the point spread function of Chandra. The best-fit 
3-/3 model is shown with red solid line, and the three components are indicated with green dashed lines. 
For the outermost /3 component, parameters of /3 and R c are fixed to the ROSAT result by Mulchaey & 
Zabludoff (1998). Intensity for the sum of the soft Galactic and the hard emission is indicated by horizontal 
magenta dashed line. The blue thick line represents the projected best-fit 2-dimensional double (3 model 
obtained with the Chandra image for the central region (r < 1'). In the bottom panel, the residuals of the 
fit are shown in unit of a. 

central region around r < V. We then varied the parameters of the inner (3 component. The 
soft background and the hard component are included as fixed constants in the fit. The fit was 
still unacceptable because of a large discrepancy in the central region. Then, we fitted with a 
3-/? model by adding a narrow component with its (3 parameter constrained to have the same 
value with the middle component, and the fit was much improved. The outermost component 
was fixed to the ROSAT value, because the IGM emission of HCG 62 is much extended over 
the XMM-Newton field of view, and we could not constrain the parameters. 

The fit results are summarized in table 3. The derived core radius, R c = 0.'48^3> for the 
middle component are consistent with the value, R c = 0/56toicL by Mulchaey & Zabludoff (1998) 
with ROSAT PSPC, while their (3 = 0.79±8;j$ is steeper than ours, (3 = 0.65 ± 0.02. This is due 
certainly to the narrow (3 model component we have introduced, which is required because 
of the superior angular resolution of Chandra ACIS-S3 than ROSAT PSPC. The obtained x 2 
value is still large, however addition of yet another (4th) (3 component does not significantly 
improve the fit, and the derived (3 values for the inner two components are close to 3.0 which 
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is obviously too steep. Therefore, we conclude that it is difficult to improve the fit better than 
the 3-/3 model case described above, due mainly to the complicated spatial structure in the 
central region. 

We next carried out a 2-dimensional fitting for the Chandra X-ray image. The ACIS-S3 
image shown in figure 3 (a) suggests that the brightest position may be slightly offset from 
the center of the extended group gas. To examine this, we fitted the 2-dimensional image in 
0.5 4 keV within r < V around HCG 62a with a 2-/3 model (narrower and wider), whose centers 
were varied as free parameters. Due to the limited photon count in the image bin, we chose the 
maximum likelihood method assuming the Poisson statistics in the fitting. The "blank-sky" 
background were included in the model as a constant surface brightness. The data and the 
best fit model are projected around the common center (HCG 62a) as shown in figure 4. The 
model (blue thick line) well describes the radial profile within r < 1', and gives an acceptable 
fit with x 2 / dof = 132.07/122. The fit parameters are summarized in table 4. We find that the 
centers of the narrower and wider components are different by 5.44" (1.6 kpc). The center of 
the narrower component is closer to the optical center (2" offset) than the wider component 
(5" offset). 

4. X-ray Cavity 

In this section, we look into the X-ray structure and the spectral characteristics of 
the cavities which are recognized in the Chandra image in figure 3 (a). We broadly consider 
two possibilities for the origin of the cavities as follows; (1) depression of the X-ray flux is 
resulted by absorption due to some intervening material, and (2) X-ray emitted gas is deficient, 
possibly being expelled by some process. We will examine the case (1) in §4.2 and (2) in §4.3, 
respectively. 

l^.t. Flux depression & temperature map 

To quantify surface brightness depression in the cavity regions, we refer to the result of 
the 2-/3 model fit of the 2-dimensional image in § 3.3. The relative deviation of the brightness 
from the 2-/3 model is color coded and plotted in figure 5. Both of the cavity regions show a 
brightness drop by as much as 50-70% at the bottom as compared with the level of surrounding 
regions. If we approximate the cavity shape by a sphere, their radii are ^ 10" (3 kpc) as indicated 
by white circles in figure 5, with the north one slightly larger. This approximation is similar 
to the one performed in B04. The distances of the cavity centers from HCG 62a are ~ 30" 
(9 kpc) for the north and ^ 20" (6 kpc) for the south one, respectively. If we take the distances 
from the central position of the wider component of the 2-/3 model, the cavities are almost 
symmetrical in their distances ~ 25" (7.4 kpc). The positions and the distances of both cavities 
are summarized in table 5. 

Next, we look into temperature distribution using hardness ratios (HR). After the point 
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Fig. 5. (a) Color-coded map of the relative deviation of the flux from the best-fit 2- (3 model based on 
the 2-dimensional fit (table 4) in 0.5-4 keV within a radius of 1' around HCG 62a with ACIS-S3 (J2000). 
X-ray contours shown by think black lines are the same in figure 3 (a). White and black circles with a 
radius of 10" represent regions where energy spectra are studied. White ones are the cavities, and black 
ones are non-cavity regions examined for comparison. The distances from HCG 62a to the centers of the 
east and west light-blue circles are 20" and 30", the same for the north and south cavities, respectively, 
(b) Residual of the observed count to the 2-/3 model along the green rectangular region (20" x 120") in 
(a), (c) Ratio of the residual counts divided by the 2- (3 model. See text §4.1 for red lines. 

source and background subtractions, we produced images in two energy bands, 0.5-0.95 keV 
(S) and 0.95-4 keV (iJ), which give nearly the same counts in both energy bands. The HR 
is defined as the ratio of the counts between these bands as HR = (H — S)/(H + S). We 
divided the region into small cells whose sizes are determined to contain at least 50 counts 
in the 0.5-4.0 keV (H + S) band, and then the HR values were calculated. The HR values 
were then converted to temperature based on the response matrices at the center of the field 
of view, assuming an absorption of TVh = 3.0 x 10 20 cm -2 and a metal abundance of 0.5 solar. 
Typical 90% confidence error on HR is ~ 0.17, which corresponds to ~ 0.2 keV above 1 keV 
and ~ 0.1 keV below 1 keV. The resultant temperature map is shown in figure 6. A clear 
temperature drop down to ~ 0.7 keV is seen at the central region, which is surrounded by a 
hotter gas with kT ~ 1.4 keV at a radius greater than V (18 kpc). The central cool region 
shows an irregular shape, with an elongation to the directions of two cavities, although the 
temperature structure does not simply correlate with the cavity regions. 
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Fig. 6. Color coded temperature map (J2000) based on the hardness ratio, HR = (H — S)/(H + £), with 
ACIS-S3, where S and correspond to the 0.5-0.95 keV and 0.95-4 keV counts, respectively. The red 
circle has a radius of V around HCG 62a, and the small circles (10" in radius) are the same as those in 
figure 5. Typical 90% confidence error on temperature is ~ 0.2 keV above 1 keV and ~ 0.1 keV below 
1 keV. 

4-2. Spectral comparison between cavity & non- cavity regions 

In order to examine the possibility that absorption may be responsible for the cavities, we 
look into the energy spectra in this subsection. We extracted energy spectra for the two cavity 
regions from the ACIS-S3 data, using circular regions with the same radius of 10" (3 kpc) as 
shown in figure 5. To compare with them, we also extracted spectra for two non-cavity regions 
with both 10" radii in the east and west at the same distance from HCG 62a to the respective 
cavities. Their precise locations are (12 h 53 m 06f87, -9°11'49"4), (12 h 53 m 05f54, -9°12'35"1), 
(12 h 53 m 07?75, -9°12'21"1) and (12 h 53 m 04?33, -9°12'20"1) (J2000) for the north, south cavities 
and east, west non-cavity, respectively. The non-cavity regions are indicated by black circles in 
figure 5, and their properties are summarized in table 5. First of all, we calculated the hardness, 
HR, as defined in the previous subsection, in these four regions. The derived values with 90% 
confidence errors were —0.130 ±0.066 vs. —0.209 ±0.052 for the north cavity vs. east non- 
cavity pair, and —0.099 ±0.058 vs. —0.1 17 ±0.046 for the south cavity vs. west non-cavity pair, 
respectively. There is no significant difference in the hardness for both pairs. However, both of 
the cavity regions suggest slightly larger hardness than the respective non-cavity regions, which 
is in the same sense as expected if absorption is responsible for the X-ray cavities. 

To be more quantitative, we conducted a spectral fit for each spectrum. All the spectra 
were binned to contain at least 30 counts in a bin. The spectra were then fitted with a two 
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Table 5. Observed and calculated properties of the north and south cavities and the east and west non-cavities. 



north cavity east non-cavity south cavity west non-cavity 

Distance from HCG 62a * 31.9" 31.9" 19.7" 19.8" 

Distance from group core + 26.5" 29.0" 24.1" 24.7" 

HR^ -0.13 ±0.07 -0.21 ±0.05 -0.10 ±0.06 -0.12 ±0.05 

Integration radius § 15.7" <— 12.6" <— 

7V c or7V nc H 2225 ±48 2857 ±58 2476 ±63 3316 ±75 

Fproj ~ Sphere or F proj ** 50.63 64.93 48.31 63.46 



* Position of HCG 62a is assumed to be coincide with the center of the narrower component in table 4. 

i" Position of the group core is assumed to be coincide with the center of the wider component in table 4. 

* Hardness ratio, HR = (H-S)/(H + S), within r < 10", where S (or H) corresponds to 0.5-0.95 (0.95-4) keV counts. 
§ Integration radius is chosen to become N c /N nc = (-Fproj — Sphere V-^proj f° r each cavity and non-cavity pair. 

II N c for north and south cavities, and N nc for non-cavities. Each value is corrected for background and exposure. 
^ ^proj — ^sphere f° r north and south cavities, and F pro j for non-cavities in arbitrary unit. 

Table 6. Results of spectral fits for the north and south cavity regions and the east and west non-cavity regions within 
r < 10" in 0.5-3 keV with a two temperature vMekal model. See text §4.2 for details. 

,excess kT\ Normi Norrri2 -^0.5— 4 keV X /dof 

(10 20 cm" 2 ) (keV) (10- 19 cm" 5 ) (lO" 14 ergs cm" 2 s" 1 ) 



phabs (TVh = 3.0 x 10 20 cm" 2 fixed) x 2-T (kT 2 = 1.4 keV & abundance fixed) vMekal model 



north cavity 

east non-cavity 

south cavity 

west non-cavity 




74+ - 04 

U -' 4 -0.05 
u -'°-0.05 

76+ - 05 

u - '°-0.06 


1 0+0.3 

L '°-0.2 

9 r+0.4 
Z '°-0.4 

1 6 +0 - 3 

2 6 +0 - 4 
z -°-0.4 


7 +0 - 4 
u - ' -0.4 

6 +0 - 5 

1 1+O.6 
J-- 1 -0.4 


4.0 ±0.3 
6.8 ±0.4 
5.2±0.3 
8.0 ±0.4 


18/14 
22/16 
18/20 
39/32 


phabs (TVh = 3.0 x 10 20 cm" 2 fixed) 


x zphabs x 2-T (kT 2 = 1.4 keV & abundance fixed) vMekal model 


north cavity 

east non- cavity 

south cavity 

west non-cavity 


18.3±|;1 

11.7l 3 2 6 3 
<2.8 


70+ - 06 

U - ' U -0.02 

T 

u - ' °-0.04 

T 


9 c+0.6 
z -°-0.4 

T 

9 r+0.4 

Z '°-0.3 

T 


7+0.3 

U -'-0.3 

T 

1 9+0.4 

T 


4.2 ±0.3 

°-°-0.3 

5.5 ±0.3 

7.6 ±0.4 


} 42/31 
} 76/53 



* Normalization for the Mekal model, Norm= J n e nndV/(4;7r(l + z) 2 D^J, where Da is the angular distance to the source. 

temperature vMekal model (2-T model in § 5) in the 0.5-3 keV band with an absorption (phabs 
in XSPEC) fixed to the Galactic value, 7V H = 3.0 x 10 20 cm -2 . Because of the limited statistics of 
those spectra, temperature of the hot component, fcT 2 , was fixed to 1.4 keV and the abundance 
of each element was fixed to the best-fit values for the 1-T model in the annulus range of 
0.4 < r < 0.6 / in table 7. See § 5 for details of the model. The fit results are summarized in the 
first four rows of table 6. 
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Fig. 7. Chandra ACIS-S3 X-ray spectra for north (black) and south (red) cavity regions, and east (green) 
and west (blue) non-cavity regions, respectively. Spectra of the north cavity and east non-cavity pair and 
the south cavity and east non-cavity pair are simultaneously fitted with common continuum parameters 
with individual excess absorption. The spectra and the best-fit models are scaled by 2, 4, 6 for south, 
east, and west regions, respectively, for clarity. The bottom four panels show the residuals of the fit. 

All of the four spectra are well fitted by the 2-T model with acceptable x 2 values at the 
90% confidence limit. The obtained temperatures exhibited quite similar values of kT\ = 0.74- 
0.76 keV, and we could not find any significant difference between the cavity and non-cavity 
spectra except for the normalizations of the cool component, Normi. It is interesting that the 
normalizations of the hot component, Norm^, are not different between them. As described in 
the next section, the hot component is probably absent from the central region of the group 
(r < 0.6'), so that the hot component seen in these four regions are mostly due to foreground 
or background emission. It is therefore suggested that the cavities lie close to the group core 
in our line of sight, and that only the cool component is pushed away from the cavities. The 
signature of hardening of HR above is due probably to this effect. As a remark, the fits are 
not acceptable when fitted with a one temperature vMekal model (1-T model), although the 
obtained temperatures are also very similar among the four spectra. 

Here, we evaluate how much A H , e xcess is needed to reproduce the observed cavity feature 
by a simple absorption at the source redshift. For this purpose, we simultaneously fitted 
the cavity and non-cavity spectra of each pair, by setting all the parameters of the emission 
spectrum to be the same including the normalization. Thus, two new parameters, A H , e xcess for 
each pair, and three new constraint on Normi, Norm 2) and kT\ are added, so that the degree of 
freedom (dof) increases by one. The pulse-height spectra and fit residuals for the four regions 
are shown in figure 7. The A H , e xcess for the individual regions are summarized in the bottom 
four rows of table 6. The obtained x 2 /dof were 42/31 for the north cavity and east non-cavity 



west X6 
east X4 

south cavity X2 
north cavity X1 




1 1 1 1 — | 
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pair, and 76/53 for the south cavity and west non-cavity pair. Therefore both fits were not 
acceptable at the 90% confidence limit, although improvement of \ 2 is marginal for the former 
pair. The derived iV H ,excess i n cavities has to be ~ 3 times larger than those in non-cavity regions 
for the necessary flux reduction. The required mass of neutral hydrogen for north and south 
cavities amount to 1.0 x 1O 9 M and 0.4 x 1O 9 M , respectively. These values are comparable to 
the observed upper limit of the Hi mass, ~ 1O 9 M , of the whole group including all member 
galaxies by Verdes-Montenegro et al. (2001) and Stevens et al. (2004). It is therefore difficult 
to attribute the cavities to the absorption by neutral gas. 

4-3. Hollow sphere model 

Since the absorption model is found to be unlikely in the previous subsection, we then 
consider an extreme case that a spherical cavity contains no X-ray emitting gas at all, namely, 
the hollow sphere model. It is assumed that the spherical cavities exist at the same distance 
from us to the group core. Even with this extreme assumption, there should be a certain 
flux observed in the projected circular cavity region because of the foreground and background 
group emission. This will give us a constraint on the physical size of the cavities. For an overall 
X-ray structure, we assume that the 3-dimensional X-ray emissivity obeys the 3-/3 model shown 
in table 3. Based on this model, we can estimate the emission from an arbitrary volume in the 
IGM and the emission from a column integrated along a certain line of sight. Here a single 
cavity is considered. The observed projected counts corrected for background and exposure 
at the circular cavity region is denoted by AT C , and the count in the corresponding non-cavity 
region by 7V nc , respectively. The non-cavity region has the same integration radius and distance 
from HCG 62a as the cavity region. Then, based on the 3-/3 model, we can calculate the flux, 
^sphere, from the cavity volume which is assumed to be filled with hot gas. The model also 
gives us the projected flux, F pro j, from the line-of-sight column passing through the non-cavity 
region. These numbers allow us to estimate the counts purely emitted from the spherical cavity 
as, N c — AT nc (Fproj — F sp here)/-Pprojj which should be zero on the assumption above. 

Based on the observed cavity flux, we can derive the upper and lower limits for the 
radius of the hollow spheres. Clearly, too small radius will give too little depression in the 
X-ray flux at the cavity region, and vise versa. We applied the observed intensity and position 
of the north and south cavities. The allowed ranges for the radii of the hollow spheres were 
computed as r = 15.7 ±0.9" and r = 12.6 ±0.8" for the north and south cavities, respectively. 
These sizes are almost consistent with the projected image as seen in figure 5 (a). The observed 
and the calculated values of AT C , 7V nc , F pro j — F sp h e re, and F pro j are summarized in table 5. The 
expected profile of the count rate deficiency with this assumption is overlaid in figure 5 (c) as a 
solid red line. The dashed red lines correspond to the radii at the ±90% confidence errors. As 
seen in this plot, the spherical assumption of cavities gives an asymmetric shape with deeper 
deficiency at the near side and shallower at the far side, when projected to the 2-dimensional 
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image. Because the observed shapes of the cavities are almost symmetrical, it is suggested that 
there must be an asymmetry in the shape and/or the density inside the cavities. Namely, the 
far-side of the cavities should be larger in size and/or weaker in the X-ray emissivity. It is also 
notable that even with this extreme assumption of hollow spheres at the same distance to the 
group core, the calculated deepness of the deficiency is nearly the same level or even slightly 
shallower than the observed ones. This fact implies that the shape of cavities is probably 
elongated in the direction of our line of sight. Another possibility is that the group gas itself 
has a flatter shape in the depth, although it is implausible considering the redshift distribution 
of the member galaxies (Zabludoff & Mulchaey 2000) and that the projected X-ray image is 
quite symmetric. 

5. Radial Profiles of Temperature & Abundance 

5.1. Consistency between Chandra & XMM-Newton 

In order to obtain radial profiles of temperature and abundance of the IGM, we examined 
energy spectra for each instrument (ACIS-S3, MOS1, MOS2, and pn) from several circular 
annuls from the center of HCG 62a. The center was taken from the narrower component in our 
2-dimensional fit carried out in § 3.3. First of all, we checked the consistency between Chandra 
and XMM-Newton using the spectra within r < 2 / around HCG 62a. Each spectrum was binned 
to contain at least 30 counts per bin to be tested with the \ 2 fit- The background data were 
taken from the same region in the "blank-sky" data for each instrument. The hard emission 
and the soft background components are also considered in the fitted model (§2.3). As for 
the XMM-Newton data, MOS1 and MOS2 spectra were summed up and simultaneously fitted 
with the pn spectrum. The energy range around the Al-K a line (1.4867 keV), the instrumental 
background, was ignored for both MOS and pn. 

We have fitted the spectra with a single temperature or two temperature vMekal model 
including the Galactic photoelectric absorption. Abundances of C, N, Na, and Al were fixed to 
be 1 solar. We grouped several elements and constrained them to have a common abundance. 
The first group contains O and Ne, the second group S, Ar and Ca, and the third group is 
for Fe and Ni. Among the other elements, abundances of Mg and Si were determined sepa- 
rately. In the case of the two temperature fit, we used the sum of two vMekal models in which 
two components were constrained to have common abundances. The actual model formula 
are phabsx(vMekal+powerlaw)+Mekal (1-T model) and phabsx (vMekal+vMekal+powerlaw) 
+Mekal (2-T model). Here, phabs represents photoelectric absorption and A^ H is fixed at the 
Galactic value of 3.0 x 10 20 cm -2 . The Mekal component is the soft background described in 
§2.3, with kT and Z both fixed at 0.3 keV and 1 solar, respectively, and with the fixed normal- 
ization. The power-law represents the hard emission described in §2.3 with its photon index, 
T = 1.5, and normalization both fixed. 
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The r < 2 / spectra are shown in figure 8 (a), and the fit results with 2-T model are 
summarized in table 7. The 1-T model was rejected at high significance with x 2 /dof= 2162/425 
(combined) for this region. Three kinds of the fit results are listed in table 7, Chandra (ACIS- 
S3) only, XMM-Newton (MOS1+2, pn) only, and the combined fit of all the instruments. We 
have confirmed that most of the results from different instruments were consistent at the 95% 
confidence limit. The abundances with ACIS-S3 were slightly larger than, but mostly consistent 
with those with XMM-Newton, while only Si gave inconsistent values at the 95% confidence 
limit. This is due probably to the systematic effect by the intrinsic Si line as seen in the 
"blank-sky" data for both Chandra and XMM-Newton. Since the Si line is most prominent in 
the IGM spectrum, the good statistics result in formally inconsistent abundance values. In the 
annular range of 2-4', spectral fit for individual instruments (ACIS-S3, MOS1+2, and pn) gave 
consistent results at the 90% confidence limit, due mainly to the poorer statistics. Based on 
these results, data from all the instruments were simultaneously fitted in the following analysis. 

5.2. Mekal vs APEC model & notes on abundance 

We also summarize results when the spectra are fitted with phabs x (vAPEC + vAPEC + 
powerlaw) + Mekal (2-T vAPEC model) in table 7. The vAPEC model (vl.3.1) gives smaller 
X 2 value especially for the Chandra spectrum, and indicates slightly higher temperatures by 
AT < 0.05 keV than the vMekal model. This is due probably to the better modeling of the Fe-L 
line complex in the vAPEC model, however it is confirmed that all the element abundances 
are consistent within 90% confidence errors between the two models. We therefore adopt the 
vMekal model in the following analysis considering comparison with previous literatures. 

Table 8 represents a result when all the element abundances in the vMekal or vAPEC 
model are determined separately for the combined fit. The errors of the Ne, Al, Ar, Ca, and Ni 
abundances are significantly larger than those of other elements, while other parameters agree 
with the values in table 7 within 90% confidence errors. Thus abundances in table 7, which are 
linked together for (O, Ne), (Mg, Al), (S, Ar, Ca), and (Fe, Ni), are essentially representing O, 
Mg, S, and Fe abundances, respectively. There is an indication that the Ni abundance is higher 
than the Fe abundance in both vMekal and vAPEC model fits. This might be explained by 
the fact that the fitted spectra are not deprojected so that they contain emission from outer 
region in the group which exhibits higher temperature than the group core as seen in § 5.4. On 
the other hand, it is known that Ni is synthesized more abundantly than Fe in SN la when 
compared in the solar unit, which will be discussed in § 7. 

5.3. Deprojection analysis 

We conducted a deprojection analysis using energy spectra from several circular annuls 
from the center of HCG 62a. The inner and outer radii of the annuls were r = 0-0. 2', 0.2- 
0.4', 0.4-0. 6', 0.6-1', 1-2', 2-4', 4-8', and 8-14', in which r represents the projected radius. The 
Chandra ACIS-S3 data were used in the range of r < 4', and the XMM-Newton data were used in 
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Fig. 8. (a) Chandra ACIS-S3 (green), MOS1+2 (black), and pn (red) spectra within r < 2' around 
HCG 62a in the 0.5-4 keV band. The three spectra are simultaneously fitted with 2-T model without 
deprojection, and the best-fit models are indicated by solid lines. Each component of the two temperature 
vMekal model for MOS1+2 is indicated by blue lines. The bottom three panels show the residuals of the 
fit. (b) Same as (a), but in the 2-4' annulus fitted with the deprojected 2-T model. The outer contributions 
and fitting models for MOS1-+-2 are indicated by blue lines. 



Table 7. Results of the projected spectral fit for r < 2' region with two temperature vMekal and vAPEC models. 
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Normalization for the Mekal model, Norm= J n e nndV / (4tt (1 + z) 2 D^J, where Da is the angular distance to the source. 

2' < r < 14'. Therefore, only the annulus of 2-4' was simultaneously examined by both satellites. 
The spectra for 0.2' < r < 1' include two cavities. Since the temperature and abundance of the 
two cavities are consistent with that in the non-cavity region as shown previously (§4.2), we 
did not exclude the cavity region and simply analyzed the annular spectra. Extraction of each 
spectrum and the treatment of the background are conducted in the same way as described in 
§5.1. 

Spherical symmetry was assumed in the deprojection procedure. Starting from the 
outermost region, we fitted the projected annular spectrum with 1-T or 2-T model. Spectrum 
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Table 8. Same as table 7, besides the element abundances are determined separately for the combined fit. 
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Table 9. Fractional contribution of the outer shells to the inner shells in the deprojection analysis. 



4-8' 2-4' 1-2' 0.6-1' 0-0.6' 

8-14' 0.373 0.075 0.018 0.004 0.002 

4-8' 0.273 0.056 0.012 0.007 

2-4' 0.273 0.048 0.026 

1-2' 0.238 0.109 

0.6-1' .... — 0.531 



of the neighboring inner region was then fitted with a model which contained contribution 
from the outer regions with fixed model parameters at their best-fit values. We repeated this 
procedure until the innermost region was reached. Because the outer contributions are fixed 
in the inner model fitting, the propagation of errors is ignored in this method. However, the 
inner regions are always brighter than outer regions, hence this effect is negligible. Table 9 
summarizes the fractional contribution of the outer shells to the inner shells. For example, the 
spectrum of the 4-8' annulus was fitted with adding 37.3% flux of the best-fit model in the 
8-14' annulus. 

Without the deprojection, the metal abundance of the group center (r < 1') is underesti- 
mated by about 30%, while the temperature does not change significantly. The sample spectra 
in the 2-4' annulus fitted with the deprojected 2-T model are shown in figure 8 (b), and the 
fit results are summarized in table 10. Figures 9 & 10 show the result of the deprojection 
analysis, and the 1-T results are quite similar to the previous ASCA result (Finoguenov & 
Ponman 1999). 

5.4- Temperature profile 

The temperature profiles obtained by both the 1-T and 2-T models with the deprojec- 
tion analysis are shown in figure 9. For the radius range 0.6 < r < 2.0', the 2-T model gives 
significantly lower x 2 values (x 2 = 78,68 for 0.6-1.0', 1.0-2.0', respectively) than the 1-T model 
(x 2 = 159, 126). The significance was tested with the F-test, and the 2-T model was preferred 
with a significance higher than 5a. Within r < 0.6', the 2-T model is preferred for the whole 
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Table 10. Results of the deprojected spectral fit with 1-T and 2-T vMekal models. The Chandra data (ACIS-S3) are used 
for 0-4', and the XMM-Newton (MOS1+2, pn) data are used for 2-14'. 
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XMM-Newton results 



2.0-4.0 1.43lg;S£ < - 67 ^-o°47 0.68lg;^ 0.54l°;;[° 0.4ll°;Jf ^-i.l 241/199 
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XMM-Newton results 



2.0-4.0 0.90lJ:£ l-83l 'J 2 0.16lJ;^ 2.13+lf 6 l-57±5:fg 1.071^ 1.02t£:°f 0.3±S:| 3.2±1 : ° 235/197 

4.0-8.0 0.24l * 1.401^1 0-13l° i 1.0412;?? 0.39l o :^ < 0.22 0.32l ^ 3.2± 3 {% 10A±]; 7 7 153/154 

8.0-14.0 0.431°:^ 0.891°:^ 0.191°$ 0.07j£J? 0.10l°:i 0.3ll°:|? O.Oel^ 23.91^-j 22.9li; : * 185/178 

combined results 

2.0-4.0 0.8918:% mlg^ 0.13±g:g 1.231^° 1.68l°£ 1.23l ; 8 70 1.05±g : |g 0.4±g : g 3.3lg : ? 402/312 

* Normalization for the vMekal model, Norm = J n e nndV/ (4tt(1 + z) 2 D^), where Da is the angular distance to the source. 

t The S abundance was fixed to the best fit in calculating errors due to an unresolved problem on XSPEC for this particular fit. 

0.0-0.6' data by Ax 2 = 26 (99.9% confidence with F-test), while the 1-T model is also accept- 
able when we split the region into smaller annuli, 0.0-0. 2', 0. 2-0.4', and 0.4-0.6'. Furthermore, 
we have also tried the 2-T model for the spectra in the range 0-0.2' and 0.2-0.4', but the nor- 
malization of the second component turned out to be nil. This is due probably to the poorer 
statistics for the smaller annuli and also to the fact that the cooler component is dominant in 
the central region. In the outer regions of r > 2', the x 2 values for the 2-T model were also lower 
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Fig. 9. (a) Temperature profiles with the 1-T vMekal models based on the deprojected spectral fit. 
ACIS-S3 (black) is used for r < 4', and MOS1, MOS2, and pn (gray) are used for 2' < r < 14'. (b) Same 
as (a), but with the 2-T vMekal models. 

than the 1-T case, but both of the fit were acceptable. For the outermost range 8-14', the 2-T 
model gave similar temperatures for the hotter and cooler components within errors, therefore 
the 1-T fit turned out to be just enough. This is because the hot component dominates the 
surface brightness. 

For the 1-T model, the temperature rises from kT = 0.7 keV at the center to kT — 1.4 keV 
at larger radii, which is consistent with the previous deprojection analysis for the ASCA data 
(Finoguenov & Ponman 1999) and the ROSAT temperature profile by Buote (2000a). This 
profile in r < 2' is also seen in the temperature map calculated from the HR (figure 6). The 
temperature with the 1-T model drops again at r > 8 / from the intermediate level of kT = 
1.4 keV, as shown in figure 9. The temperature drop in our data is much steeper than those 
reported by Finoguenov k Ponman (1999) with ASCA and by Buote (2000a) with ROSAT. 
The ROSAT temperature at r > 8' is kT = 0.9-1.0 keV, while ours is kT = 0.66 ± 0.03 keV. The 
influence of the background has been tested by changing the background level within ±5%, and 
we found that the best fit temperature varied in 0.64-0.70 keV. We have looked into the effects 
of point source contamination and the influence of hard and soft components, though none of 
them gives significant change in temperature. 

For the 2-T model, temperature of the cooler component is nearly constant at ~ 0.7 keV, 
which is very close to the central temperature obtained with the 1-T fit. Temperature of the 
hotter component is also nearly constant at ^ 1.4 keV up to r ^ 8 / . These features are consistent 
with the ASCA results, which indicate two temperature components at 0.7 keV and 1.4 keV 
based on projected spectra within r < 3 / (Buote 2000b). The cool component is concentrated 
in the center and dominant within 1', while the hot component extends out to ^ 4 / . The extent 
of the cool component will be examined in § 6 more quantatively. 
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Fig. 10. Abundance profiles of O (with Ne) (a), Mg (b), Si (c), and Fe (with Ni) (d) using the 1-T 
(crosses) and 2-T (diamonds) vMekal models with the deprojection analysis. The ACIS-S3 (black) is 
used for r < 4', and MOS1, MOS2, and pn (gray) for 2' < r < 14'. The solar photospheric value of 
[Fe/H] = 4.7 x 10 -5 (Anders and Grevesse 1989) is adopted for the Fe abundance. The solid and dashed 
lines correspond to the best fit relations of the deprojected O and Fe profile of the 2-T results, respectively. 

5. 5. A bundance profile 

The metal abundances were derived from the deprojected spectra based on the 1-T and 
2-T fits. The Fe, Si, and Mg abundances are around one solar or more in r < 4 / On the 
other hand, the O abundance is always lower than the solar value (table 10 and figure 10). 
We could poorly constrain the S abundance due to the limited statistics. The 2-T fit gave 
significantly higher abundance than the 1-T fit in the intermediate range of 1.0 / < r < 4.0', 
while similar abundances between 1-T and 2-T are obtained for the inner (r < 1.0') and outer 
(4.0' < r) regions. Such discrepancy is often seen in the spectral fit, because the abundance 
takes the lowest value at the temperature where the line emissivity is the highest. This situation 
corresponds to the case of 1-T fit, and the 2-T model tends to shift the temperature away from 
such a high emissivity position. The ASCA spectra of HGG 62 are fitted with the 1-T Mekal 
model by Finoguenov & Ponman (1999), and the present results for Fe, Si, and Mg with 1-T 
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fits are in good agreement with the ASCA deprojection results within the 90% confidence limit 
at r > 1', 

In figure 10, the solid and dashed lines correspond to the best fit relations of the de- 
projected O (a) and Fe (d) profile of the 2-T results, respectively. A simple linear model of 
Z = —ar + b is assumed for each profile. These best fit relations are overlaid to all the four plots 
for comparison. O abundance of ~ 0.3 solar is significantly lower than the other elements. It 
also indicates a flatter profile, though significance is low due to the low abundance of oxygen. 
The decline of Fe, Si, and Mg abundances with radius previously reported by Finoguenov & 
Ponman (1999) with ASCA, is also confirmed with Chandra and XMM-Newton. Our result 
shows that Fe, Si, and Mg abundances are 1-2 solar at the center and drop to about 0.1 solar 
at r ~ 10 / . These abundance levels are about twice larger than the values derived with the 
ASCA spectra (~ 0.6 solar by Finoguenov & Ponman 1999). This is due primarily to the much 
better angular resolution of the Chandra X-ray telescope (0.5") than ASCA (3 7 ) at the very 
central region. In addition, our 2-T model gives higher abundances in the intermediate range 
of 1.0' < r < 4.0'. As described in §5.4, the 2-T model is significantly needed in this region, 
so that the ASCA abundances are likely to be underestimated. We also note that many au- 
thors recently take the solar abundance to be the value given by Grevesse & Sauval (1998), i.e. 
[Fe/H] = 3.2 x 10 -5 , which is obtained from the measurements of solar system meteorites. The 
solar photospheric value of [Fe/H] = 4.7 x 10 -5 (Anders and Grevesse 1989), which is adopted 
in our analysis, gives the Fe abundance approximately by a factor of 1.47 smaller. Considering 
this effect, all the Fe, Si, and Mg abundances have similar value around 1.2 solar at the center 
(r <0.6'). 

6. Mass Profiles 

6.1. Formulation of two-phase gas 

Based on the radial distributions of the cool and hot components, we will derive mass 
distributions of gas and dark matter. Since the estimated gas mass contained in a single cavity 
volume, if it is filled up, is ~ 6% of the mass in the shell of 7" < r < 32", the assumption of 
the spherical symmetry gives relatively small errors. However, because two temperatures are 
required for each shell within the deprojected radius of 0.6 / ^4.0 / as seen in the previous section, 
we need an additional assumption between the hot and cool components. Here we assume a 
pressure balance between the two phase, i.e. P gas = n\kT\ = n 2 /cT 2 , where ni, n 2 , T\ and T 2 
are the cool and hot gas densities and temperatures, respectively, at each 3-dimensional radius 
of R. Such a pressure balance is previously adopted by Ikebe et al. (1999) for the Centaurus 
cluster, and by Xue et al. (2004) for the RGH 80 galaxy group. 

We also introduce a volume filling factor, f{R) ) of the cool component, namely a frac- 
tional volume, Vi = fV^ in the total volume, V is filled with the cool gas, and V2 = (1 — /) V 
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is filled with the hot gas. This means that the cool and hot gas are not completely mixed, 
instead the cool component is somewhat patchy or localized with rather an irregular shape. 
Such a hypothesis may be supported by the 2-dimensional temperature map (figure 6) and the 
existence of the cavities. 

Because the normalizations of the 2-T vMekal models are expressed as Normi = C\ 2 n\V\ 
and Norm 2 = C\ 2 n\V 2 using a certain common constant, C12, the volume filling factor, 
can be solved under the pressure balance as, 



f(R) = 



1 + 



T 2 \ 2 Norm 2 ^ ' 



(2) 



vTi/ Norm 1 

We have calculated this formula for each shell in 0.6-1.0', 1 .0-2.0', 2.0-4.0', and 4.0-8.0' using 
the best-fit values in table 10 to plot f(R) against the 3- dimensional radius R in figure 11 (a). 
This plot clearly indicates that the cool component is dominant at the central region, while it 
occupies only < 1% at R > 4'. We have fitted these four points with a 3-dimensional /5-model 
function, f(R) = [1 + (R/ R cJ ) 2 }~ 3l3f/2 , and obtained R cJ = 0.43' ±° H and f3 f = 0.60l^. 

Considering the projection effect, we can calculate the volume occupied by the cool gas 
at each 2-dimensional ring of 0i n -6 OVLt in unit of radian, as 

v x = / \ tt Di (e 3 out - C) (i - # r 2 ) 3/ 7 (i - o*), (3) 

where 9 r = 0[ n /9 out , and Da = 61 Mpc is the angular diameter distance to HCG 62. Because 
the vMekal normalization is defined as Norm = / n e nndV/ (47r(l + z) 2 D]^) ) the electron density 
of the cool gas, n e i, is computed as 



n el = ^1.2x4n(l + zfDl/V u (4) 

assuming n e — 1.2 nn for a fully ionized gas with hydrogen and helium mass fraction of X = 0.7 
and Y = 0.28. The electron density of the hot gas, n e2 , is similarly derived by replacing / into 
(1 — /). The ion density including helium is n\ = 0.92 n e , therefore the gas pressure is calculated 

as 

P gas = 1.92 n el m = 1.92 n e2 kT 2 . (5) 

The derived gas pressure P gas , cool or hot gas temperature kT\ or kT 2) and cool or hot 
electron density n e i or n e2 are plotted in figures 11 (b)-(d). The 2-T fit results are adopted 
for points in the range of 0.6-4' and the 1-T fit results for others. The systematic errors when 
the background level is increased or decreased by ±5% are considered in the error bars for the 
outermost two points. It is supposed that blue points belongs to the cool component and red 
points to the hot one. 

We must be careful in dealing with these plots because data points are not inde- 
pendent with each other. We first fitted the gas pressure plot by combining both the cool 
and hot data points, with a 3-dimensional double /3-model, P gas = Pi [1 + (i?/i? c ,Pi) 2 ] _3 ^ Pl ^ 2 + 
P 2 [1 + (R/ i? c ,p 2 ) 2 ] _3 ^ P2 ^ 2 - Due to the small number of data points, we have fixed i? c ,Pi = 0.1' and 
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(3p 1 =0.65 using the best-fit values in table 3 obtained with the radial surface brightness profile 
fit. The fitted parameters are P 1 = 58 ±130 eV cm- 3 , P 2 = 18 ±32 eV cm" 3 , i? c?P2 = 0.60 ± 1.44', 
and f3p 2 = 0.38 ±0.23. The best-fit model is indicated by a dashed line in figure 11 (c). The 
green curves represent each component of the double /3-model. 

We then fitted the cool gas temperature with a power-law model, kTi(R) = aR b , 
and the hot gas electron density with a 3-dimensional /3-model and a constant, n e2 (R) = 
n e 2,o [l + (R/Rc,n e2 ) 2 r 3Pne2/2 + C n e2 . The best-fit values are a = 0.80 ± 0. 17 keV, 6 = 0.09±0.20, 
n e2 ,o = 0.007 ±0.020 cm" 3 , i? c?ne2 = 1.4 ±7.2', (3 ne2 = 1.2 ±5.8, and C Ue2 = (0.49 ± 0.48) x 
10 -3 cm -3 . The best-fit models for kT\ and n e2 are drawn by dashed lines. The model func- 
tions for kT 2 and n e i are derived using the relation kT 2 = P gas /^e2 and n e i = P g ^/{kT 2 )^ which 
are indicated by another dashed line in each panel. Although errors of these parameters are 
quite large, these functions can reproduce the observed properties of the cool and hot gases 
fairly well. 

6.2. Dark matter and gas mass 

The gas mass density p gas is expressed as 

p gas = 1.92 /imp [fn el + (1 - /)n e2 ], (6) 

where \i = 0.62 is the mean molecular weight, and m p is the proton mass. Assuming the 
hydrostatic equilibrium, the total integrated gravitational mass, M <R) within the 3-dimensional 
radius of R is given by 

R? dP 

<R PgasG dR ' ^ ^ 

in which G is the gravitational constant. The differential mass density, M(i2), is given by 

, w 1 dM <R , x 

Figures 11 (e) and (f) show the integrated and differential mass profiles (black lines) as a 
function of the 3-dimensional radius in unit of arcmin or kpc. We also overlaid the gas mass 
(gray) for the hot (red) and cool (blue) components in the same panels, by integrating the gas 
density, p gas . 

However, in the two-phase model under the pressure balance, the hot gas is lighter in 
mass density than the cool gas, and becomes buoyant. The hot gas would escape from the group 
core if it is not enclosed by other mechanisms (see §8 in detail). In such a case, hydrostatic 
equilibrium may have been broken, hence we have also calculated the integrated gravitational 
mass by treating the hot and cool components separately. Namely, 

M h<R = ~^^^ pi = 1.92/xm p n el , 
M 2 , <i? = -^^, p 2 = 1.92/xm p n e2 , 
M <R = fM l9<R + (l-f)M 29<R . (9) 
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Fig. 11. (a) Volume filling factor f(R) of the cool component, (b) cool or hot gas temperature kT\ or 
kT2, and (c) gas pressure P gas , (d) cool or hot electron density n e i or n e 2, plotted against the 3-dimensional 
radius R. Blue points represent the cool component and red points the hot. In panels (b)-(d), 2-T fit 
results are adopted in the range of 0.6-4' and 1-T fit results for others. The systematic errors when the 
background level is increased or decreased by ±5% are considered in the error bars for the outermost two 
points. The dashed lines in panels (a)-(d) and green curves in panel (b) represent the best fit models. See 
text in details, (e) A black line indicates the total integrated gravitational mass, M<#, estimated from the 
best fit models. The solid line is calculated by eq. (7), and dashed line is by eq. (9). Red and blue lines 
represent the gas mass of the hot and cool components, respectively. The O and Fe mass contained in 
either of the gas is plotted by cyan and pink lines. T^g stellar mass estimated from the R band photometry 
of the HCG 62a galaxy by Tran et al. (2001) is indicated by a green line. See text for details, (f) Same as 
(e), but for the differential mass density profiles. 



The gravitational mass indicated by dashed lines in figures 11 (e) and (f) are calculated using 
eq. (9). There are only small differences between the two at the group center; the latter is 
larger by about 15% at R = 0.1' (1.8 kpc). 

We have encountered one severe problem in these plots. The integrated total mass 
saturates at R ~ 6 / (100 kpc), and the differential mass becomes even smaller than the gas 
mass. This is physically unrealistic, suggesting that some presumption(s) we have supposed 
might not be realized around this radius. The direct source of the problem is caused by the 
fact that the hot gas temperature, fcT 2 , drops very steeply from ~ 1.5 keV to 0.64 keV in this 
radial range. Therefore, the electron density, n e , almost saturates, while the gas pressure, P gas , 
decreases monotonously. This might in part be affected by the smaller field of view in our 
observations than the extended group emission. We cannot deproject the outermost shell in 
8-14', so that X-ray flux in this shell is slightly overestimated. Nevertheless, the major origin 
of the problem is due to the steep temperature drop. It should be examined more precisely 
by the Suzaku satellite, which has lower and more stable background than XMM-Newton, and 
much superior low-energy sensitivity than ASCA. 

Another possibility, which is astrophysically more interesting, is that the hydrostatic 
equilibrium is not reached there. The flattening of the electron density suggests an outflow of 
the hot gas. The steep temperature drop can be reversely recognized as the temperature rise in 
the boundary region between the cool and hot gas. This suggests a shock heating of accreting 
ICM, and/or remnants of past lifted cavities (see §8). The apparent inconsistency in the total 
mass may therefore indicate traces of the dynamical evolution. 

Aside from the above problem, the derived gravitational mass is about M <R = 5 x 1O 12 M 
within 100 kpc, close to the value given by Pildis et al. (1995), M tot = 2.9 x 1O 13 M within 
r < 15 / (270 kpc) with the ROSAT PSPC observation assuming a single-phase gas model. We 
also plot the O (cyan) and Fe (pink) mass contained in the whole gas. The simple linear fit 
plotted in figure 10 is adopted to calculate the O and Fe mass. The derived gas and Fe masses of 
M g8iS <R = 1.6 x 1O 12 M and M Fe<R = 4 x 1O 8 M within r < 300 kpc are consistent with the result 
from ASCA and ROSAT by Finoguenov k Ponman (1999) of M gas<R = (1.7 ±0.1) x 1O 12 M 
and M Fe<R = (3.8 ±2.0) x 1O 8 M . However, our values are based on the extrapolation out to 
the XMM-Newton field of view, so that they have larger errors by a factor of ^ 2. 

6.3. Stellar mass 

The stellar mass in figures 11 (e) and (f) indicated by green lines is estimated from 
the R band photometry of the HCG 62a galaxy by Tran et al. (2001), assuming the spherical 
symmetry and the surface brightness profile obeying the de Vaucouleurs r 1 / 4 law, 

E^oclO" 3 - 33071 ^ 174 - 1 ], (10) 

where r e ff = 32" is the effective radius, or the half light radius. Because the r 1 / 4 law is for the 
surface brightness, deprojection for the 3-dimensional radius must be performed. We utilized a 
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numerical table computed by Young (1976). The mass-to-light ratio, M stdbT /L B = 8 (M /Lb >0 ), 
is assumed, and we use B — R = 2.0 mag by Hickson et al. (1989), and ^ = 0.224 and t4r = 0.139 
after NASA/IPAC Extragalactic Database (NED) to calculate the B band absolute magnitude, 
Lb = L R + [B — R) — A B + A R = —19.46 mag. Using the value L B ^ = 5.48 mag, the stellar 
mass of the HCG 62a galaxy is calculated as M star = 7.6 x 1O 1O M . 

The B magnitude within the effective diameter of 59.90" isophotes (table 1) is also given 
by Hickson et al. (1989), however, we found that M star calculated from this value was by about 
a factor of four larger than the value calculated above, if we assume that r e ff is the same for both 
bands. Then, the stellar mass would significantly exceed the total gravitational mass at the 
central region. We therefore have adopted the stellar mass based on the R band photometry. 
Tran et al. (2001) also report that there is 19% of disk component which can be fitted by 
E(r) ocexp(— r/r&) for HCG 62a. This gives asymmetric image residual flux Ra = 19%, and the 
total residual fraction of light Ft? = 18%. Therefore, uncertainty on the stellar mass estimation 
is at least ^ 20%. 

We notice several features in the mass profiles. The gas mass of the cool component is 
taken over by the hot component at ~ V (18 kpc). The stellar mass is overcome by the gas 
mass at ^ 2.5' (45 kpc), indicating that the cool component is very concentrated in the group 
center. The region of high metal abundance almost coincides with the volume dominated by 
the stellar mass, which supports the natural view that stars are responsible for the production 
of excess iron and silicon around the group center. 

7. ICM properties 

The temperature structure in HCG 62 can be characterized by a mixture of hot (~ 1.4 
keV) and cool (^0.7 keV) components. As shown in figure 11, the cool component is dominant 
within r < 0.8 / and then a cool and hot mixture appears in r = 0.8-4'. The hot component 
is dominant in the outer region. The two-phase nature seems to be preferred than the single- 
phase structure such as that seen in M 87 (Matsushita et al. 2002). However, this apparent 
two-phase does not necessarily mean that the cool and hot gas co-exist at the intermediate 
region, instead it probably represents that they are patchy and/or have irregular shapes. As 
seen in figure 6, the cool region is elongated from northeast to southwest, roughly corresponding 
to the direction of the two cavities. In this sense, HCG 62 group is different from NGC 1399, 
NGC 5044 and RGH 80 galaxy groups (Buote 2002; Buote et al. 2003a; Xue et al. 2004), in 
which strong evidence for a multi-phase gas is suggested at the central region. On the other 
hand, T coo i is similar to these groups, and T coo i c± is recognized among them. It is claimed 

that the temperatures of the cool component seen in these three groups are close to the kinetic 
temperature of the stars, although the stellar velocity dispersion is not known for the HCG 62a 
galaxy. 

Mass density of the cool component exceeds the hot-component density within about 
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10 kpc, in which the stellar mass density is orders of magnitude higher. This again implies 
that the cool component is probably connected with the stars concentrated around the central 
galaxy HCG 62a. We also note that the cool component mass density is always smaller than the 
stellar mass density. On the other hand, the hot component exceeds the stellar mass density 
around r ~ 2 / (30 kpc), indicating that the gravitational potential of the galaxy group is traced 
by the hot component. Its temperature (1.4 keV) is also typical for groups of galaxies. 

In the outer region (r > 8'), the temperature of the hot component drops fairly sharply. 
The radius corresponds to about 0.15r vir , with r vir « 1.1 (T/1.5 keV) 1 / 2 Mpc. This causes a 
peak in the temperature profile at r « 5 / (90 kpc). Similar temperature profiles have been 
observed in other galaxy groups such as RGH 80, NGC 2563, and NGC 5044 (Xue et al. 2004, 
Mushotzky et al. 2003, Buote et al. 2003a). The drop of temperature suggests either that the 
dark matter is confined within this small radius or that the gas is yet to be heated in this region. 
In either case, there is a certain boundary of the group around 0.2r v i r , and it is suggested that 
these groups are young and forming systems. 

The heavy elements are enriched by SN la and SN II (Tsujimoto et al. 1995). The 
former process dominantly yields Fe group (Iwamoto et al. 1999), and the latter produces lighter 
elements like O, Ne and Mg (Thielemann et al. 1996; Nomoto et al. 2006). Both supernova 
contribute to Si and S. Since O is produced mainly by SN II, the flatter O abundance profile 
compared with those of Si and Fe implies enhanced contribution of SN la in the central region 
and/or that the shallow potential of HCG 62 is unable to confine SN II products which should 
have been supplied to the intracluster space in the form of galactic winds. The products of 
SN la, on the other hand, are considered to be brought in to ICM by gas stripping. 

The ratio, Mg/O, is 3.3 ±2.2 times solar within V (18 kpc) based on the 2-T fit. Though 
the error is large, this value is similar or somewhat larger than those in other groups: 2.5 ±0.4 
solar for NGC 5044 (Xu et al. 2002), 1.3 ±0.2 for NGC 4636 (Tamura et al. 2003), both 
measured with RGS, and ~ 2 solar in RGH 80 (Xue et al. 2004). In M 87, this ratio is 1.3±0.1 
solar within 20 kpc (Matsushita et al. 2003). Since both O and Mg are mainly produced by 
SN II, the marginally high Mg/O value may reflect the difference in the stellar initial mass 
function (IMF). Theoretical calculations (Nomoto et al. 2006; Thielemann et al. 1996) predict 
that a lighter-mass progenitor of SN II synthesizes enhanced Mg compared with O, therefore 
HCG 62 might have had the stellar IMF with smaller number of massive stars (M > 2OM ) 
than our galaxy and other groups above. However, the dependence of the Mg/O ratio on the 
progenitor mass is typically ~ 20% at most, so that we have to think about other possibilities 
for the Mg/O ratio greater than ~ 1.2. 

If excess Fe in the center is caused by SN la activity of HCG 62a, it should be more 
extended at least than the extent of the central cool component which is the gas directly 
bound by the central galaxy. The possible enhancement of the Ni/Fe ratio at r < 2 / as seen in 
§5.1 supports this scenario. For example, the CDD2 model in Iwamoto et al. (1999) predicts 
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the Ni/Fe ratio of 1.8 solar. The region of high metal abundance almost covers the volume 
dominated by the stellar mass, which supports the natural view that stars are responsible for 
the production of excess Fe and Si around the group center. It seems to range even to the 3D 
radius of ~ 6' (100 kpc) in figure 10, where the hot gas mass exceeds the stellar mass as seen in 
figure 11 (f). This suggests either that the Fe production has occurred in a wider region than 
the present location of the galaxy or that there have been an outflow of Fe. The presence of 
cavities in HCG 62 suggests that part of the central metal-rich gas may have been lifted from 
the galaxy. 

The observed abundance gradients indicate that there have been no strong mixing oc- 
curred in the core region. We also note that there is no significant change of temperature or 
abundance across the cavity regions. These features imply that the process of cavity creation 
causes mild, subsonic gas motion occurring in a fairly limited volume. This implication is con- 
sistent with the result by Briiggen & Kaiser (2002) who showed that mixing by buoyant bubbles 
gave relatively weak impact on the metallicity gradients based on numerical simulations. 

8. Supporting mechanisms of X-ray cavities 

B04 have systematically studied 16 clusters, 1 group (HCG 62), and 1 galaxy (M84), in 
which prominent X-ray surface brightness depressions (cavities or bubbles) are observed. They 
find that a mechanical (kinetic) luminosity seems to correlate with the 1.4 GHz synchrotron 
luminosity. However, its ratio ranges widely between a few and several thousand, and they have 
concluded that the radio luminosity is an unreliable gauge of the mechanical power of the AGN 
jets. Dunn & Fabian (2004) and D05 have studied 21 clusters, 3 galaxies, in which HCG 62 
is not included, and find that the ratio of an energy factor, /C, to a volume filling factor, T, 
shows a large scatter, 1 < JC/J 7 ^ 1000, for active cavities associated with radio lobes, and that 
it becomes even larger for ghost cavities. The factor, /C, accounts for the additional energy 
from relativistic particles accompanying the electrons (e.g., protons). The parameter of T 
represents the volume filling factor of the relativistic particles, and is supposed not to vary 
far from unity. Typical value of JC in literature is JC = 100, according to the measurements of 
cosmic-rays around the solar system. However, there is no direct evidence indicating that such 
a high energy density is really carried by protons. The reason why JC/ J 7 varies so largely from 
cluster to cluster is a mystery, which might imply that there are several ways in supporting 
mechanism and/or formation of the X-ray cavities. 

In this section, we examine whether this non-thermal pressure support scenario is re- 
alistic or not for HCG 62. In §8.1, we summarize our X-ray results, as well as optical and 
radio observations. In §8.2-8.4, standard indices of cavities are calculated for HCG 62, and 
difficulties in the non-thermal support scenario is considered in §8.5. In §8.6 and §8.7, other 
possibilities of the supporting mechanism and formation of cavities are investigated. 
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8.1. Characteristics in X-ray, optical, and radio 

We have confirmed the X-ray cavities reported by Vrtilek et al. (2001,2002) in the north- 
east and southwest regions of HCG 62 in the Chandra image as shown in figures 3 (a) and 5. 
As described in § 4.2, absorption is unlikely to be the origin of cavities. There are no significant 
spectral differences between cavity and non-cavity regions, nor any traces of the shock-heated 
gas around the cavities, as seen in the temperature map (figure 6). According to our hollow 
sphere model analysis (§4.3), both the cavities should be aligned side by side with the group 
core to a fairly good degree, namely, the projected distances to the group core are close to their 
real ones. It is also suggested that the shape of cavities is probably elongated in the direction 
of our line of sight, and that the side of the cavities which is farther from the core should be 
larger in size and/or weaker in the X-ray emissivity. 

With regard to the central HCG 62a galaxy, we found no evidence for the AGN activity 
in the X-ray data. The best-fit spectrum of HCG 62a is a thermal emission rather than a 
power-law one. No point sources were recognized at the center of HCG 62a in both soft and 
hard bands. The upper limit of the AGN emission can be placed around L x ^ 10 39 erg s _1 
(0.5-4 keV). We also found that the location of HCG 62a was slightly shifted from the group 
core by about 5" (1.5 kpc) on the projected sky image (table 4). This positional shift may 
imply the effect of gas stripping and/or the dynamical motion of the galaxy. 

The ASCA has detected a spatially extended (~ 10') excess hard X-ray emission above 
~ 4 keV in the HCG 62 group (Fukazawa et al. 2001), which is supposed to be due to the 
relativistic electrons with Lorentz factor 7 ~ 10 3 -10 4 and/or sub-relativistic particles. This 
may also have some relation to the X-ray cavity, however, we could not confirm nor reject this 
result due to the higher non X-ray background of Chandra and XMM-Newton than ASCA. 

In the optical band, Coziol et al. (1998) and Shimada et al. (2000) have detected weak 
[Nil] and [OI] emission lines in the spectrum of HCG 62a, and classified it as a low luminosity 
AGN (LLAGN). On the other hand, Coziol et al. (2004) have assigned the lowest activity 
index of —5 (quiescent, intermediate and old stellar populations) to HCG 62a, according to 
the equivalent width measurement of the Ha absorption lines. They have also distinguished 
HCG 62 group among 27 compact groups of galaxies as type C, which comprises groups with 
high velocity dispersions and are dominated by elliptical galaxies with no activity, presumably 
corresponding to the later stage of the evolution. 

In the radio band, as shown in figure 3 (a), weak radio emission at 1.4 GHz is detected 
around HCG 62a, however its luminosity is as small as I/ radio = 1.8 x 10 38 erg s _1 (10 MHz-5 
GHz), smallest in table 1 of B04. The radio emission does not show a clear association with the 
cavities, though the angular resolution of 45" FWHM is not sufficient to see detailed structures. 
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8.2. Gas pressure & non-thermal pressure 

We take the distances to both cavities from the group core to be i? cav = 25" (7.4 kpc) 
and their radii to be r cav = 13. 5" (4.0 kpc) for simplicity. This approximation can be justified 
by the similarity between north and south cavities measured from the group core (table 5) with 
the center position determined from the center of the wider component in 2-dimensional 2-/3 
fit (table 4). Assuming those parameters, the ICM gas pressure at i? cav is calculated to be 
^Pgas = 17 eV cm -3 , from eq. (5). 

The non-thermal pressure of relativistic particles can be estimated from the radio in- 
tensity. Here, we assume that half of the radio flux comes from a single cavity, namely, 
^radio = 9 x 10 37 erg s" 1 (10 MHz-5 GHz) according to B04. Assuming that the radio emission 
is due to synchrotron radiation by relativistic electrons, the total non-thermal pressure can be 
calculated as P t ot = lCP e + J 7 Pb, where P e is a pressure of the relativistic electrons, and Pb is 
a pressure of the magnetic field, following the convention adopted by D05. 

The relativistic electron pressure is calculated to be P e = C e L ra di £> -3 / 2 /V (e.g., Govoni 
& Fertti 2005), where B is a magnetic flux density, V = §7rr 3 av is the volume of the cav- 
ity, and C e is a constant depending on the spectral index a of the radio emission, as 



C e = ^27ra5c 9 /(27re 7 ) (^~ 1/2 - z/ 2 " 1/2 )/(logz/ 2 - log^) = 1.4 x 10 8 [cgs] for a = 1, v x = 10 MHz, 
and ^2 = 5 GHz. Therefore, 



JCP e = 3.3 x 



^radio 



. , f-g-V 

9xl0 37 cgs/ ^4kpc/ 



.100, 

On the other hand, the magnetic pressure is 



eV cm 3 . 



Therefore, the total non-thermal pressure, Ft = fcPe + J 7 Pb, takes the minimum at a certain 
value of B = £> eq , the so-called equipartition condition, which is calculated as 




^T,eq — 7^e,eq — ~ ft,eq — 5.8 X 

4 _ 12 

\ 7 / m \ 7 

eV cm 3 

9xl0 37 cgsy \4kpcy 

Thus the derived pressure at the equipartition condition is less than half of the ICM 
gas pressure, while the required equipartition magnetic flux density is typical for radio lobes 
of radio-loud AGNs (Kataoka & Stawarz 2005). They have also found that the equipartition 
condition is achieved at least for the radio lobes of 40 radio galaxies. The discrepancy between 
P T ,e q and P gas may be explained by the underestimation of JC/J 7 ^ 100 or L radio ^9 x 10 37 erg s _1 , 
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or by non-equilibrium situation. To make a balance between them, JC/J r = 690 is needed. This 
value is consistent with those obtained by D05 for active bubbles in other clusters. It is also 
plausible that the radio intensity is getting dimmer in time due to the synchrotron cooling of 
relativist ic electrons. We also note that the I/ radio value used for the calculation is probably 
overestimated, because the observed radio intensity is likely to include emission from the core 
regions. Radio observation with higher angular resolution is desired. 

8.3. Time scales 

For discussion of cooling and non-equilibrium effects, comparison of several time scales 
is important. First of all, synchrotron electron (7 e > 10 4 ) must be long-lived since the cavity is 
not fueled now from the central AGN. The synchrotron cooling time is calculated as 

H 3 c 5 or / b V 2 fie Y 1 „ 

Since rims of the cavities have not been shock heated, the cavities are supposed to have expanded 



at a velocity less than the sound speed, v s = y 7/cT/ (/xm p ) = 425 km s _1 , where we have taken 
kT = 0.7 keV, 7 = 5/3 and fi = 0.62. Therefore, age of the cavities must be longer than the 
expansion time of 

which is shorter than i sync . If we assume JC/J r = 690, i sync can be shortened by a factor 
of (/C/JT/100) -4 / 7 = 0.33, and becomes comparable to t s . The time scale for the cavities to 
collapse when internal pressure has disappeared is also considered to be about t s . Such hollow 
cavities are buoyant and rise up outwards, even when pressure balance between inside and 
outside of the cavities is conserved. Churazov et al. (2001) have given the terminal velocity as 
where S = 7rr^ av is the cross section of the bubble, g = GM <Rcav / R 2 ^ is the 
gravity at the bubble, and C = 0.75 is the drag coefficient. The travel time of the cavities to 
the current position is estimated to be 



^buoy Ft* 




-Kcav \ / ^cav \ / M <R C&V \ 



Myr, 



V7.4kpcy V4kpcy V 2xl ° llM ©/ 
which is comparable to or somewhat longer than t sync = 8 Myr with JC/J 7 = 690. The refill time 
of the cavity is also given by McNamara et al. (2000) as, 




Refill = 2R C3W x /r cav /(GM <Rcav ) = 31 x 

' M <Rc&v 




V7.4kpcy V2xlO n M 
These considerations indicate that the cavity age of t s or tbuoy ~ 10 Myr is much shorter than 
the group age of ~ Gyr. Our time scales are consistent with B04, but in which they use R cscv 
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instead of r cav . 
8.4- Energetics 

On the assumption that past AGN activity has produced the two cavities in about 
(tbuoy/2), we can estimate the required mechanical (kinetic) power of the AGN jets. A work to 
generate two cavities, V^ mec h = 2P gas y, divided by (ibuoy/2) is called the mechanical luminosity 
(B04), and calculated to be 



It is supposed that this level of the AGN activity must continue for about — t huoy <t< —0.5t huoy , 
and that it became inactive since —0.5 ibuoy < t. This value is comparable to the X-ray emission 
typically observed for LLAGN, although our upper limit on the X-ray luminosity is much 
lower, L x < 10 39 erg s" 1 (0.5-4 keV). The observed radio emission around HCG 62a, L radio = 
1.8 x 10 38 erg cm -1 , is also much lower than Z/ mec h, as pointed out by B04. 

Because there are no evidence for the strong AGN activity at present time, the central 
AGN, if exists, should have made a final outburst of total energy ~ 10 57 erg within ~ 20 Myr 
ago. Clearly, a single supernovae cannot account for this size of energy. Such an absence of 
strong X-ray or radio emission at the core is also noticed in NGC 4636 (I/ ra dio = 1.4 x 10 38 erg s _1 ; 
L x < 2.7 x 10 38 erg s" 1 ; Ohto at al. 2003; O'Sullivan et al. 2005), which shows a disturbed X- 
ray halo containing cavities associated with small-size jets. This might suggest that the AGN 
activity continues only for a short time scale, which is difficult to understand within the popular 
paradigm of AGN with a steady accretion disk. Formation of ghost cavities with an impulsively 
episodic activity of <C 10 5 yr is discussed by Wang & Hu (2005). Considering captures of red 
giant stars by a super massive black hole (M B h > 2 x 1O 8 M ), its feedback energy can amount 
to 2.4 x 10 52 erg with a frequency of a few 10 -5 yr -1 , which is too small to supply sufficient 
energy for the formation of cavities. It is probably true that HCG 62a also contains a massive 
black hole of M B h ~ 1O 8 M , because observations of the centers of nearby early- type galaxies 
show almost all have massive black hole (Tremaine et al. 2002). However, we have to consider 
switching (on — > off) of the AGN activity. 

8.5. Difficulties in non-thermal pressure support by jets 

In § 8.2-8.4, we have examined the properties of the cavities and the central AGN from 
the point of view that the non-thermal pressure of relativistic particles were provided through 
symmetrical jets of the past AGN activity. However there are several difficulties in this scenario. 

As seen in §8.2, K/J 7 = 690 is needed to balance the ambient gas pressure with the 
internal non-thermal pressure at the equipartition condition. This value is in proportion to the 
inverse of L ra dio, which is probably overestimated because significant fraction of Z/ ra dio should be 



£mech = 4P gas y/t buoy = 1.8 X 10 42 X 
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attributed to the core region instead of the radio robes (figure 3 (a)). Then JC/J 7 = 690 becomes 
still larger, which is probably unrealistic. One possibility is that the relativistic electrons are 
fading due to the synchrotron loss, which may be justified by the fact that i sync under the 
equipartition magnetic field of B eq = 17 /iG (JC/J r = 690) is comparable to the estimated cavity 
age of tbuqy- The relativistic electrons have much shorter synchrotron cooling time than that of 
protons, hence apparently large JC/ ' T might have been attained. In this case, the cavities are 
starting to collapse and losing the internal pressure support. 

It is notable that both the two cavities observed in HCG 62 are pretty circular when 
projected in the sky (figure 3 (a)). The relative deviation image of figure 5 is remarkably 
smooth around the group core, and there are no obvious trails toward the cavities. This fact 
is quite difficult to understand considering the scenario that a pair of radio robes produced by 
symmetrical jets from the AGN has pushed away the IGM. In practice, clusters or galaxies 
hosting cavities usually show irregular or filamentary structures around the cluster core and 
cavities. The most prominent example is those of M 87 (Churazov et al. 2001; Young et 
al. 2002), in which the "trails" of the rising two radio bubbles are clearly seen in the X-ray 
image. Furthermore, our hollow sphere analysis suggests that the shape of cavities is probably 
elongated in the direction of our line of sight. It appears to be difficult for the symmetrical 
jets to make an elongation in such direction. We also note that shock heating is unlikely to be 
the origin of the cavities, because there is no significant evidence of heating at the edge of the 
cavities. 

These considerations indicate that there might exist another supporting mechanism 
and/or formation scenario of the X-ray cavities, at least for HCG 62. In the following sub- 
sections, we consider these possibilities. 

8.6. Another supporting mechanism — hot gas clump — 

One possible supporting mechanism is a clump of hotter gas than the surrounding ICM 
(kT = 0.7 keV). For example, if we assume that the temperatures in the cavities are higher than 
that of ICM by three times, namely Th c = 3T, the required density, n e ^ c , is three times smaller. 
The emitted X-ray is roughly in proportion to n 2 ehc \/Th~c = 0.2nlVT, hence the observed X-ray 
intensity can be as low as 20%. Since we could have detected such hot emission if its flux were 
^ 20% of the F S phere in table 5, we can derive the lower limit of the hot clump temperature, as 
T hc >3T. 

Schmidt et al. (2002) have done a similar discussion and ruled out volume-filling X-ray 
gas with temperature below 11 keV for a cavity in the Perseus cluster. With regard to the 
origin of the hot gas clump, it is suggested that the jets may intrinsically contain protons, or 
that they may have captured ambient thermal protons possibly shock-heated in the very initial 
phase of the cavity formation (D05). If we observe in higher energy band, such hot emission can 
be detected as a hard tail of the spectrum. The hard T = 1.5 power-law component observed 
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with ASCA (Fukazawa et al. 2001) might originate from this kind of hot thermal emission, 
although the detected sky area of the power-law component is much more extended than the 
cavities. 

If such a hot clump really exists, it undergoes a cooling by thermal conduction. Assuming 
kT hc = 3kT~2 keV and n e ^ c = n e /3 ~ 10" 3 cm -3 , the thermal conductivity for a non-magnetized 
plasma is given by Spitzer (1962) as, 

,, _ in 3Q ( kT hc y ( n e , hc \ AnA\ 2 

where In A is the Coulomb logarithm. Therefore, the cooling time of the cavities is roughly 
calculated to be 

y>2 
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which is much shorter than the estimated cavity life span of tbuoy 

It is claimed that the thermal conductivity may become 5-10 times smaller than the 
Spitzer value under the turbulent magnetic fields (e.g., Chandran & Maron 2004), although 
^cond seems to be still smaller considering this effect. On the other hand, from the observational 
point of view, we do see this kind of temperature variations in clusters. The most prominent 
example is the "cold front", first reported by Markevitch et al. (2000) for A2142. Ettori & 
Fabian (2000) have pointed out that it requires the classical Spitzer thermal conductivity to 
be reduced at least by a factor of 250-2500. Markevitch et al. (2003) also find that ~ 40 times 
reduction is needed for A754. For the magnetic field expected in the radio robe (B ~ 10 /xG), 
the electron and proton gyro radii are by 11-12 orders of magnitude smaller than their Coulomb 
mean free paths, therefore the effective conductivity strongly depends on the topology of the 
field. If the magnetic field encloses the hot clump like a cage, the hot clump may survive for 
more than ~ 20 Myr. 

8.7. Another formation scenario — galaxy motion — 

As we found in §3.3 with the 2-dimensional 2-/3 model fit, the location of HCG 62a 
(narrower component) is slightly offset from the group core (wider component) by about 5" 
(1.5 kpc) on the projected sky image. This fact inevitably leads to the idea that the HCG 62a 
galaxy is moving around the group core. It is natural to consider that HCG 62a is performing a 
pseudo-Kepler motion, gradually decreasing its distance to the group core by a dragging force. 
From this point of view, the two cavities might be a piece of the trail of the HCG 62a orbit. 
Here, we assume that the orbital plane of HCG 62a is nearly in parallel to our line of sight, with 
a circular orbit, and that its orbital radius is close to the distance, i? ca v, of the two cavities to 
the group core. 
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It is interesting that the two cavities are located almost at the same distance from the 
group core, whereas the south cavity is by 1.6 times closer to HCG 62a than the north cavity. 
As indicated by a cross in figure 3 (a), HCG 62a is shifted from the group core indicated by a 
plus mark (center of the wider component) roughly toward the direction of the south cavity. 
These two facts support the assumption above. The measured redshift of HCG 62a is consistent 
with that of the group (table 1). On this assumption, the rotation speed is calculated to be 



^rot = jGM <Rcav /i? cav = 340 km s x , which is comparable to the sound speed, v s = 425 km s 



Even though it is at least shorter than the group age, is much longer than other time 
scales, and the "tunnel" seems to be filled relatively quickly. In this scenario, the ICM gas 
would have experienced the encounter with the HCG 62a galaxy several times periodically. If 
some processes, e.g., freezing the plasma with magnetic field or the pressure support with hotter 
gas, have slowed the collapse of the cavities, this effect may be worthwhile consideration. In 
terms of the energetics, the kinetic energy of the galaxy motion amounts to ~ 10 59 erg, therefore 
it can supply sufficient energy to produce the cavities by depositing 7% of the kinetic energy 
per orbit. Such a motion of the galaxy also would have played an important role on the mixing 
and the metal enrichment of the ICM, which has been discussed in the previous section. 

9. Conclusion 

• We have carried out a detailed study on the hot-gas properties of the group of galaxies 
HCG 62, based on the data from Chandra and XMM-Newton. We confirmed the two 
cavities located almost symmetrically around the central galaxy. 

• The size of spherical hollow cavities are constrained from the surface brightness structure 
to be 12"-1 7". The agreement with the observed angular size suggests that the gas density 
in the cavity is very low, less than 1/3 and consistent with zero. 

• The spectral fit indicated that the cavities were not caused by X-ray absorption. The 
observed temperature in the cavity region is consistent with that in the surrounding region. 

• The spectrum within 4 / from the center requires two temperatures: 0.7 keV and 1.4 keV. 
The cool component is centrally concentrated, narrower than the hot component, suggest- 
ing its association with the central galaxy HCG 62a. 

• The mass profiles were obtained for the gas and stars. The hot component is much more 
extended than the stars, and thought to trace the gravitational potential of the galaxy 
group. 

• The gravitational mass density drops steeply at about 5 / from the center. This is caused 
by the observed sharp drop of the temperature. There is a possibility that these regions 
are not in the hydrostatic equilibrium. 




The period of rotation is 
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• Abundance of is ~ 0.3 solar, ~ 3 times less abundant than Fe and Si, and shows a flatter 
profile. The shallow potential of HCG 62 is unable to confine the SN II products which 
should have been escaped in the form of galactic winds. The marginally higher Mg/O 
ratio of 3.3 ±2.2 implies steeper IMF. 

• Abundances of Fe and Si show concentration in the central region, and a high Ni/Fe ratio 
is suggested. These results are consistent with that they are synthesized by SN la in the 
central galaxy. 

• The non-thermal energy density necessary to support the cavity implies JC/J r = 690, namely 
almost 700 times larger energy than that of electrons needs to be contained the cavity. 
The lack of the central AGN or the trailing radio feature seems to suggests that the origin 
of the cavity in HCG 62 may not be directly related to AGN activities. 

• We looked into alternative scenarios for the cavity creation. A clump of very hot gas and 
fast motion of the central galaxy were considered, but more observational evidences are 
necessary to perform a quantitative evaluation. 
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